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摘 要 : 天 体 随时 间 的 光度 变化 能 够 反映 其 物理 性 质 ,高 精度 测 光 技 术 的 研究 帮助 我 们 更 好 地 分 析 与 
探究 天 体 的 光 变 情况 ， 对 天 体 的 演化 研究 有 重要 意义 。 本 文 介绍 了 相关 的 测 光 工 具 与 技术 ,探索 了 
把 MaxIm DL 软件 与 photutils 测 光 库 相 结合 的 测 光 技 术 及 其 应 用 ,提出 了 连续 图 像 中 各 星 像 测 光 数 
据 快 速 匹 配 的 解决 方案 。 通 过 分 析 中 国 科 学 院 云南 天 文 台 Im 光学 望远镜 拍摄 的 M35 星团 图 像 ， 本 
文 发 现 了 相对 测 光 中 存在 精度 回报 (precision premium) 现 象 ， 即 相近 的 两 颗 星 进行 相对 测 光 具有 更 高 
的 测量 精度 ， 且 该 现象 的 显著 程度 与 大 气质 量 相关 。 针 对 于 相对 测 光 的 精度 回报 现象 ， 本 文 提出 了 


高 精度 相对 测 光 的 参考 方法 。 
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天 体 的 光度 测量 工作 主要 研究 天 体 的 光度 随时 间 的 变化 情况 ， 光 度 变化 能 够 反映 出 天 体 的 物理 
性 质 ， 对 探究 天 体 的 演化 规律 有 重要 意义 。 例 如 ， 文 献 [从 小 行星 的 光 变 周期 中 分 析出 自转 周期 和 
形状 ， 文献 [2] 研 究 了 柯 伊 伯 带 天 体 的 光 变 周期 ,分 析 了 其 反照 率 等 物理 属性 ， 文 献 [3] 利 用 中 长 周期 
光 变 研究 类 星体 的 物理 属性 。 然 而 ， 人 们 在 天 体 的 高 精度 光度 测量 中 面临 许多 挑战 。 由 于 地 球 大 气 
影响， 天 体 的 光度 经 过 地 球 大 气 后 存在 一 定 的 消光 和 折射 效应 ， 对 天 体 的 测量 精度 有 着 较 大 的 影 
响 。 由 于 相近 (天 空 平面 内 角 距 离 较 小 ) 的 恒星 之 间 相 对 位 置 变 化 较 小 且 通 过 的 大 气相 似 , 我 们 通常 使 
用 相对 测 光 的 方法 进行 测量 。 然 而 , 太阳系 天 体 (特别 是 近 地 小 行星 ) 相 对 于 恒星 位 置 变 化 较 快 且 经 过 
不 同 的 大 气质 量 , 绝对 测 光 (全 天 测 光 ) 则 是 一 种 更 好 的 选择 。 在 地 面 望远镜 观测 中 ， 人 们 可 以 通过 该 
方法 准确 地 求 出 大 气 消光 、 夜 间 零 点 等 系数 ， 以 降低 大 气 的 影响 。Landolt 提供 了 标准 星 的 详细 参数 
与 信息 所 ， 为 绝对 测 光 研究 商定 了 良好 的 基础 。 尽 管 如 此 ， 绝 对 测 光 的 精度 不 如 相对 测 光 。 为 了 消 
除 地 球 大 气 对 光度 测量 的 影响 ， 随 着 科技 的 发 展 和 进步 ， 人 们 把 探测 器 发 送 到 空间 中 (例如 哈 勃 空间 
望远镜 和 Gaia 卫星 )， 测 光 精 度 有 了 显著 的 提高 。 

本 文 主要 研究 基于 地 面 观测 的 光度 测量 技术 和 精度 回报 现象 。 在 地 面 观测 中 ， 为 了 提高 位 置 和 
光度 的 测量 精度 ， 通 常 我 们 选择 目标 在 子午 线 附 近 观 测 ， 以 减少 大 气 造 成 的 影响 。 在 天 体 测量 
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HERKIMER", Lin 等 人 基于 M35 星团 对 该 现象 进行 统计 分 析 ， 把 位 置 测量 的 精度 回报 现 


象 进行 特征 化 2 。 在 相对 测 光 中 ， 人 们 主要 关注 参考 星光 度 的 稳定 性 而 较 少 地 关注 目标 星 与 参考 星 


在 天 空 平面 内 的 角 距离 。 若 相对 测 光 中 也 存在 精度 回报 现象 ， 即 假设 目标 星 存 在 许多 不 同 相 对 距离 
的 候选 参考 星 ， 且 这 些 参 考 星 的 光度 相近 、 稳 定 ， 则 采用 相对 距离 较 小 的 参考 星 会 有 较 高 的 测量 精 


度 。 基 于 中 国 科学 院 云南 天 文 台 lm 光学 望远镜 观测 的 M35 星团 资料 ， 本 文 将 探究 相对 测 光 的 精度 


回报 现象 。 


* 基 金 项 目 : 国家 自然 基金 (11873026，1 


1273014); 国家 自然 科学 天 文 


育 计划 项 目 资助 ， 作者 简介 : MWEE, y 


J’ 信 息 系 统 项 


KES ae a 


EE 点 项 


管理 师 ， 硕 了 


2 


(U1431227); 暨南 大 学 优秀 本 科 推 免 


E 科 研 创 新 培 


FE 研究 生 . 研究 方向 : 图 像 处 理 . Email: bifengguo@qq.com 


Chinax Vv iv 合 人 a EH 


本 文 的 第 2 部 分 介绍 了 现 有 的 测 光 工具 及 其 技术 ， 提 出 了 把 Maxim DL 软件 与 photutils 库 相 合 
的 测 光 技术 及 其 应 用 。 第 3 部 分 统计 和 分 析 了 M35 星团 在 不 同 高 度 角 中 相对 测 光 的 精度 与 相对 距离 
之 间 的 关系 ， 发 现 了 相对 测 光 中 存在 精度 回报 现象 ， 提 出 了 高 精度 相对 测 光 的 参考 方法 。 最 后 一 部 
分 进行 了 全 文 的 总 结 。 
2 ”光度 测量 工具 与 关键 技术 
2.1 光度 测量 工具 

随 着 信息 化 与 工业 化 的 发 展 ， 目 前 有 许多 开放 的 光度 测量 工具 和 软件 。 例 如 IRAF 能 够 较 好 地 
处 理 图 像 并 进行 测 光 ， 但 需要 在 Linux 环境 下 进行 配置 和 操作 。Aperture Photometry Tool (APT) 是 一 
款 孔 径 测 光 的 教学 软件 ， 能 够 方便 地 进行 孔径 测 光 、 求 取 生长 曲线 等 站 。SExtractor 被 广泛 运用 在 深 
度 巡 天 的 目标 源 提取 和 测量 中 。MaxIm DL 除了 可 以 方便 地 进行 批量 平 场 处 理 、 图 像 质 量 评 佑 、 图 
像 堆 琶 外 ， 也 能 方便 地 进行 孔径 测 光 包 。 后 来 ， 随 着 计算 机 技术 的 发 展 和 Python 编程 语言 的 兴起 
PyRAF(https://pypi.org/project/pyraf/) 和 photutils(https://photutils.readthedocs.io/en/stable/index.html) 成 
为 了 方便 测 光 的 Python 库 。 
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2.2 MaxIm DL 与 photutils 相 结合 的 测 光 技 术 


Photutils 是 基于 Python 语言 的 测 光 库 ， 也 是 Astropy 的 附属 软件 包 ， 用 户 可 以 通过 编程 方便 地 
= 调用 其 目标 检测 、 光 度 测量 等 函数 。DaoStarFinder 和 IRAFStarFinder 是 该 库 函数 中 用 于 目标 检测 的 
N 两 个 函数 ， 分 别 调用 了 DAOPHOT 和 IRAF 的 搜 星 程序 。Photutils 库 中 提供 了 孔径 测 光 和 点 扩散 函 
数 测 光 两 种 方法 。 本 文 主要 使 用 和 探究 孔径 测 光 方法 ， 即 photutils.aperture 库 。 在 该 库 中 ， 浮 点 型 数 
值 的 孔径 半径 的 设置 为 孔径 测 光 提供 了 更 高 的 精度 ， 而 Maxim DL. APT 等 测 光 软件 的 孔径 半径 数 
值 只 能 设置 为 整数 。 此 外 ，aperture_photometry 函数 中 的 exact 方法 能 够 精确 地 处 理 孔径 测 光 中 包含 
S 像素 的 分 数 部 分 ， 相 比划 分 子 像素 (subpixel) 方 法 有 着 更 高 的 测量 精度 。 

= 然而 , 使 用 Photutils 对 地 面 的 观测 资料 进行 光度 测量 的 过 程 中 也 存在 一 些 问 题 。 在 目标 检测 ( 搜 
星 ) 函 数 中 ， 无 论 我 们 使 用 DaoStarFinder 还 是 IRAFStarFinder JE, (RABI (threshold 参数 ) 和 星 像 
的 半 高 全 宽 (FWHM 参数 ) 就 是 必 选 的 输入 参数 。 通常 , 像素 闵 值 我 们 选择 几 倍 (例如 3-5 倍 ) 背 景 像素 
的 标准 差 作为 输入 ， 而 做 巡天 研究 工作 时 需要 把 该 参数 设置 得 更 小 。 要 准确 地 搜索 出 观测 资料 中 的 
大 部 分 星 像 ， 传 入 的 半 高 全 宽 参 数 需要 与 每 幅 图 像 相近 。 当 然 ， 我 们 希望 每 个 晚上 观测 资料 的 半 高 
全 宽 是 稳定 的 ， 这 也 需要 天 气 与 设备 的 配合 。 但 实际 上 ， 在 地 面 观测 中 ， 即 使 用 同一 台 望 远 镜 连续 
观测 的 图 像 ， 星 像 的 半 高 全 宽 也 有 一 定 的 波动 。 因 此 在 测 光 程序 批量 处 理 时 ， 若 使 用 一 个 大 致 的 半 
高 全 宽 作为 目标 检测 函数 的 参数 输入 ， 则 一 些 图 像 中 大 部 分 的 星 像 不 能 被 检测 出 来 。 对 于 此 类 情况 ， 
我 们 可 以 把 半 高 全 宽 与 设 定 值 相 差 较 大 的 图 像 进 行 单独 处 理 或 编写 自 适 应 调整 半 高 全 宽 参 数 的 相关 
算法 。 由 于 上 述 方法 操作 或 实现 过 程 较 复杂 , 本 文 提出 了 与 Maxim DL 软件 相 结合 的 方法 .MaxIm DL 
软件 (本 文 使 用 5.24 版 本 ) 在 图 像 堆 登 过 程 中 提供 了 批量 图 像 质 量 评估 的 功能 ， 可 以 方便 地 导出 每 幅 
图 像 的 半 高 全 宽 、 圆 度 等 信息 (如 图 1 左 子 图 所 示 )。 当 我 们 使 用 MaxIm DL 输出 每 幅 图 像 的 半 高 全 宽 
数据 作为 目标 检测 函数 参数 的 输入 时 , 可 以 较 好 地 解决 图 像 之 间 星 像 半 高 全 宽 波 动 较 大 的 搜 星 问题 。 
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我 们 在 实践 中 发 现 ， 目 标 检测 函数 中 传 入 的 半 高 全 宽 参 数 取 Maxim DL 导出 的 数值 再 加 上 0.5 有 着 


更 好 的 搜 星 效 果 。 


在 光度 测量 的 过 程 中 ， 分 析 和 统计 出 视 场 中 所 有 星 像 的 光 变 情况 对 变星 巡天 、 星 像 特 征 分 析 有 


着 重要 意义 。 如 何 把 各 图 像 的 星 像 进 行 相互 匹配 是 统计 视 场 9 
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PhP 所 有 星 像 的 光 变 的 关键 问题 。 对 于 望 


远 镜 指向 和 取向 基本 稳定 的 观测 资料 来 说 , 获取 主 帧 (一 般 选择 质量 较 好 的 图 像 ) 所 有 星 像 的 位 置 和 光 


度 信息 就 能 较 好 地 与 其 他 图 像 匹配 (例如 设置 一 定 的 像素 阔 值 和 光度 阔 值 )， 


从 而 方便 地 画 出 目标 的 光 


变 曲线 。 然 而 ， 观 测 时 望远镜 指向 或 跟踪 不 稳定 (例如 每 幅 图 之 间 平 移 5px)， 统 计 出 视 场 中 所 有 星 像 


的 光度 变化 是 不 容易 的 。 本 文 结合 Maxim DL 软件 对 该 问题 提出 了 一 种 快速 的 解决 方案 。 首 先 ， 使 


用 MaxIm DL 的 对 齐 功 能 可 以 方便 地 导出 每 幅 图 相对 于 第 一 幅 图 在 像素 坐标 的 偏 移 量 (如 图 1 右 子 图 
P 减 去 相应 的 偏 移 量 (XShift 和 YShift)。 此 时 ， 我 们 


所 示 )。 然 后 ， 在 每 一 幅 图 探测 出 来 的 星 像 坐标 


可 以 消除 由 于 望远镜 指向 不 稳定 带 来 的 视 场 偏 移 , 能 够 
但 如 果 需 要 对 各 星 像 中 的 更 多 信息 (例如 色 指 数 、 视 差 等 ) 进 行 分 析 和 统计 研究 , 我 们 需要 采用 基于 星 
表 匹 配 的 方法 ， 但 该 方法 需要 花费 更 长 的 时 间 。 


方 


更 地 进行 各 图 像 


中 星 像 的 匹配 和 统计 工作 。 


FWHM, Roundness, Intensity, Image 
5:905; 0-162; 37273:6, ‘QL flat.fi 
T0343 0.207, 38790.5, 02 flat.fi 
5.748 0.129, 37537-1; 03 flatsfi 
5. 989, 0.148, 36169.7, 04 flat.fi 
Gat Say 0.154, 35987-3; 05 flat:fi 
6.024, 0.047, 35652.1, 06 flat.fi 
6.555 0.168, 33000-9, 07 flat.fi 
6.398, 0.148, 33755.1, 08 flat.fi 
6.048, 0.106, 34309.1, 09 flat.fi 
6.016 0.064, 34107.5, 10 flat.fi 
6-365; 0.034, 31404-1; 11 flat.fi 
6.536, 0.156, 32321.8, 12 flat.fi 
6.138 0.043, 31529.0, 13 flat.fi 
6.952, 0.125, 27071.4, 14 flat.fi 
T.047, 0.073, 26810.6, 15 flat.fi 
6.115, 0.018, 26743.7, 16 flat. 
7,007, 0:083; 22335.1, - 17 flats 
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¥Shift, Theta, Image 
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ye i 360.00, 09 flat.fit 
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-3.51, 0.00, 11 flat.fit 
-2.47, 0.00, 12 flat.fit 
-3.95, 0.00, 13 flat.fit 
-4.54, 0.00, 14 flat.fit 
-6.47, 0.00, 15- flät-fit 
-7.28, 0.00, 16. flat.£it 
-6.98 0.00; -17 flat.fit 


图 1 左右 两 幅 子 图 分 别 展 示 了 使 用 MaxIm DL 软件 导出 的 图 像 质量 评估 文件 和 图 像 对 齐 文件 


Fig.1 The left and the right panel show the quality evaluation and image alignment file exported by MaxIm DL respectively 


结合 Photutils 测 光 库 和 MaxIm DL 软件 ， 


测 过 程 中 CCD 芯片 不 产生 较 大 角度 旋转 。 
3 ”相对 测 光 中 的 精度 回报 现象 


本 节 提 出 了 由 于 星 像 半 高 全 宽 变 化 较 大 时 部 分 图 像 中 


目标 不 能 被 检测 的 解决 方法 。 此 外 ， 本 节 也 提出 了 一 种 快速 的 本 地 星 像 匹 配 与 统计 方法 ， 但 要 求 观 


受 大 气 消 流 的 影响 ，CCD 芯片 中 星 像 的 位 置 与 实际 位 置 存在 不 同 的 偏差 。 由 于 视 场 内 相近 的 两 


颗 星 受 大 气 油 流 的 影响 相似 ， 星 像 的 位 置 测 


量 中 存在 精 


立 兽 测量 有 着 更 高 的 精度 。 在 相对 测 光 中 ， 我 人 


度 回 


报效 应 ， 即 星 像 之 间 角 距离 较 小 时 相对 
] 猜 测 在 同一 视 场 内 受 不 同 大 气质 量 的 影响 ， 各 星 像 


存在 不 同 程度 的 消光 。 若 相近 的 两 颗 星 受 大 气 消 光 的 程度 相似 ， 则 和 它们 也 有 较 高 的 相对 测 光 精度 。 
在 本 节 中， 我们 将 统计 M35 星团 的 相对 测 光 精度 与 相对 距离 的 关系 ， 探 讨 测 光 中 的 精度 回报 现象 。 
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3.1 观测 资料 选取 


找到 两 是 相互 靠近 且 光 度 不 变 的 天 体 是 探究 该 现象 存在 的 最 好 方法 。 然 而 ， 我 们 很 难 找到 角 距 


离 变 化 较 大 且 光 度 不 变 的 天 体 。 对 于 太阳 系 小 行星 来 说 ， 虽 
上 自转、 相位、 卫星 等 因素 的 影响 也 会 发 生 》 
一 般 的 恒星 作为 实验 的 观测 资料 ， 由 
系 是 比较 困难 的 。 此 外 ， 也 有 不 少 恒 4 
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在 3.2 节 中 详细 讨论 对 于 这 些 f 


3.2 资料 处 理 与 数据 分 析 


O 


使 用 了 


百 星 的 处 理 。 


EC 度 的 变化 ， 不 适合 月 


H, RIIK 


is 
pr 


星团 


然 它 能 够 与 恒星 相互 靠近 ， 但 小 行星 受 
日 于 探究 该 现象 的 存在 性 问题 。 若 以 
于 恒星 之 间 的 相对 角 距 离 几 乎 不 变 ， 找 到 测 光 精度 与 距离 的 关 
的 亮度 是 变化 的 。 在 实验 
资料 。 星 团 中 大 量 的 恒星 让 我 们 方便 地 统计 
该 现象 是 普遍 存在 的 ， 而 不 是 只 发 生 在 某 个 天 区 或 个 别 天 体 上 。 因 此 ， 使 用 


] 星 团 作为 实验 的 观测 
距离 和 测 光 精度 的 关系 。 此 外 ， 星 团 的 统计 能 够 反映 


E 够 较 好 地 探究 相 


E 变 化 或 受 其 他 天 体 光 度 影 


国 科学 院 云 南天 文 台 lm 光学 望远镜 (详细 参数 如 表 1 所 示 ) 括 
测 光 精度 回报 现象 的 探究 。 我 们 在 不 同 的 高 度 角 对 M35 星团 
县 体 的 观测 资料 情况 如 表 2 所 示 。 


表 1 望远镜 及 CCD 详细 说 明 
Tab.1 Specifications of the Telescope and CCD 


Item Parameter 
Focal Length 1330cm 
Diameter of Primary Mirror 101.6cm 
CCD Field of View 16’ X16’ 
Size of Pixel 15u X15u 
Size of CCD Array 4096 X 4112 
Scale ~0.234’’/px 


表 2 观测 资料 概要 


Tab.2 Observations Overview 


响 较 大 的 恒星 ， 我 们 


摄 的 M35 星团 图 像 进行 


进行 了 观测 , 观测 时 望远镜 的 指向 稳定 ， 


Date(y-m-d) Frames/Filter Altitude(”) Exptime(s) FWHM(px) 
2020-12-08 16/I 83-89 40 5.5-6.5 
2021-01-15 21/1 83-78 40 4.1-5.3 
2021-01-15 20/1 50-43 50 4.6-6.0 
2021-01-15 17/1 36-31 50 5.7-7.0 


对 于 观测 资料 的 预 处 至 


L RIEA 
选 出 半 高 全 宽 较 小 (例如 小 于 6px， 即 


图 像 帧 数 是 经 过 质量 评估 第 选 后 的 数量 。 在 实际 的 测 光 
第 选 ， 因 为 该 方法 可 能 把 关键 的 光 变 数据 剔除 。 在 本 实验 中 ， 探 究 相 对 测 光 的 精度 回报 现象 的 存在 


需要 假定 每 颗 恒 星 的 光度 者 
地 减少 由 于 大 气 的 不 稳定 、 


是 不 变 ， 
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日 MaxIm DL 4x4 


NJ 


也 需要 尽 可 能 地 把 光度 
焦距 等 因素 造成 的 影响 。 


做 平 场 校正 ， 用 其 
\ 于 1.4 角 秒 ) 且 圆 度 较 好 (例如 小 于 0.2) 
PF， 我 们 要 慎重 使 用 


中 


像 的 质量 评估 功能 簿 


该 方法 


发 生变 化 的 星 像 剔除 ， 该 方法 能 够 有 效 


的 图 像 。 表 1 中 列 出 的 


进行 图 像 质量 的 
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我 们 使 用 孔径 测 光 的 方法 进行 实验 探究 。 由 于 星团 中 恒星 的 星 像 比较 密集 ， 使 用 环形 区 域 来 估 
计 各 星 像 的 背景 值 常常 会 受到 附近 星 像 的 干扰 。 因 此 ， 本 文 在 每 幅 图 像 的 各 星 像 中 都 采用 一 个 固定 
的 天 空 背景 值 。 具 体 地 ， 我 们 先 对 每 幅 图 像 背 景 像素 吻 除 了 306 ， 然 后 采用 “3-2” 公 式 (3 倍 中 值 -2 
倍 均值 六 "进行 天 空 背景 值 的 估计 。 接 着 ， 我 们 使 用 photutils 测 光 库 中 提供 的 函数 进行 搜 星 ， 结 合 
MaxIm DL 软件 进行 各 星 像 光度 的 匹配 和 统计 ， 有 具体 见 本 文 2.2 节 。 因 为 我 们 需要 假设 每 颗 恒 星 在 每 
组 图 像 的 光度 都 是 不 变 的 ， 故 以 相对 测 光 仪器 星 等 在 每 组 图 像 的 标准 差 进 行 衡量 。 由 于 上 暗 星 测 光 精 
度 较 低 ， 我 们 选择 信 咯 比较 高 的 亮 星 来 统计 星 像 的 相对 测 光 精度 与 相对 距离 之 间 的 关系 。 这 些 亮 星 
分 布 在 视 场 不 同 的 位 置 ， 我 们 让 这 些 亮 星 两 两 进行 较 差 测 光 ， 统 计 每 组 图 像 中 星 像 相 对 测 光 精度 和 
与 参考 星相 对 像素 距离 之 间 的 关系 。 这 时 ， 若 每 组 图 像 中 一 共 选 择 出 n 颗 亮 星 ， 每 颗 被 选择 的 亮 星 
都 会 作为 其 他 n-1 颗 星 的 参考 星 ， 一 共有 n(n-1) 组 测 光 精度 与 距离 的 关系 。 
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图 2 左 子 图 展示 了 2020 年 12 A 8 日 资料 测 光 精 度 与 相对 距离 之 间 的 关系 ， 统 计 了 仪器 星 等 为 11-13 的 数据 。 相 对 
像素 距离 采用 公式 JAYAY ， 其 中 Ax Ay 为 星 像 中 心 的 像素 坐标 之 差 (下 同 )。 右 子 图 是 使 用 “标准 差 均 值 剔 除 ” 
方法 (设置 立 值 为 0.0055) 后 的 结果 ， 红 色 实 线 是 使 用 肾 函 数 拟 合 的 结果 。 


Fig.2 The left panel shows the relationship between the relative distance and photometric precision of the observations on 


December 8, 2020, and the data of instrumental magnitude range 11-13 are counted. The relative pixel distance adopts the 


formula Ax? +Ay? » where Ax and Ay are the difference of pixel coordinates of star image center (same hereinafter). The 


right panel is the result of using "standard deviation mean elimination" method (threshold set as 0.0055), and the red solid 


line is fitted by power function. 

然而 ， 以 信 噪 比 或 星 等 作为 标准 筛选 出 来 亮 星 的 光度 也 不 一 定 是 稳定 的 。 若 某 一 恒星 是 变星 、 
双星 系统 或 受 附 近 亮 星光 度 影响 较 大 时 ， 该 恒星 的 测 光 精度 也 会 比较 低 ， 也 不 适合 用 于 统计 和 分 析 
相对 测 光 的 精度 回报 现象 。 因 此 ， 我 们 需要 尽 可 能 找 出 这 些 测 量 精 度 较 低 的 星 ， 并 把 它们 吻 除 。 
们 知道 ， 若 一 颗 星 的 光度 在 一 组 图 像 中 发 生 了 较 大 的 变化 ， 它 不 管 与 哪里 参考 星 进 行 相对 测 光 ， 都 
会 有 较 低 的 精度 (标准 差 数值 较 大 )。 本 文 提出 了 “标准 差 均值 剔除 ”的 方法 来 剔除 测量 精度 较 低 的 恒 
星 。 第 一 步 ， 我 们 先 计算 出 每 颗 恒 星 与 视 场 中 其 他 恒星 的 星 等 标准 差 ， 则 每 蜂 星 都 能 找到 n-1 组 这 
样 的 标准 差 。 第 二 步 ， 求 出 每 颗 星 对 于 n-1 组 标准 差 的 均值 。 第 三 步 ， 设 定 一 个 效 值 (一 般 设 置 为 星 
像 精度 集中 的 上 方 附 近 ， 如 图 2 左 子 图 所 示 )， 把 标准 差 的 均值 高 于 该 阅 值 的 星 别 除 。 使 用 该 方法 能 
够 有 效 地 吻 除 精度 不 高 的 恒星 ， 方 便 统计 和 探究 测 光 精度 回报 现象 的 存在 。 下 图 展示 了 2020 年 12 
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月 8 日 拍摄 的 M35 星团 图 像 对 于 测 光 相对 精度 与 像素 距离 之 间 的 关系 ,其 中 采用 了 仪器 星 等 为 11-13 
的 恒星 。 左边 子 图 是 吻 除 之 前 的 结果 ,右边 子 图 采用 “标准 差 均值 剔除 ”的 方法 ,设置 阔 值 为 0.0055 


剔除 后 的 结果 ， 可 以 看 到 能 够 有 效 地 剔除 精度 较 差 的 星 。 
接着 , 我 们 使 用 同样 的 方法 处 理 了 2021 年 1 月 15 日 共 3 组 在 不 同 高 度 凶 拍摄 的 M35 星团 资料 ， 
5 所 示 。 


PA 


以 探究 大 气质 量 对 测 光 精度 回报 现象 的 影响 ， 结 果 如 图 3- 


0.012 


Mag(12-14,threshold = 0.008) 


© 
己 
© 
te] 


Precision(mag) 
Residual(mag) 


0.003 


0 1000 2000 3000 0 1000 2000 3000 


Distance(px) Distance(px) 


图 3 左 子 图 展示 了 2021 年 1 月 15 日 资料 (高 度 角 为 36~31”， 仪 器 星 等 为 12-14mag) 测 光 精 度 与 相对 距离 之 间 的 关 
系 ， 其 中 剔除 阔 值 设置 为 0.008， 右 子 图 是 使 用 窜 函 数 拟 合 后 残 差 随 相对 像素 距离 之 间 的 分 布 情况 。 


Fig.3 The left panel shows the relationship between the relative distance and photometric precision of the observations on 


j 


January 15, 2021 (the altitude is 36~31 deg and instrumental magnitude is 12-14mag), in which the elimination threshold is 


set to 0.008. The right panel shows the distribution of the residual fitted by power function. 


Mag(11-13, threshold = 0.006) 


Precision(mag) 
Residual(mag) 


0 1000 2000 3000 0 1000 2000 3000 


Distance(mag) Distance(px) 


图 4 左 子 图 展示 了 2021 年 1 月 15 日 资料 (高 度 角 为 30 ~ 43”， 仪 器 星 等 为 11-13mag) 测 光 精 度 与 相对 距离 之 间 的 关 
系 ， 其 中 剔除 阔 值 设置 为 0.006， 右 子 图 是 使 用 肾 函 数 拟 合 后 残 差 随 相 对 像素 距离 之 间 的 分 布 情况 。 


Fig.4 The left panel shows the relationship between the relative distance and photometric precision of the observations on 


sal 


January 15, 2021 (the altitude is 50 ~ 43 deg and instrumental magnitude is 11-13mag), in which the elimination threshold is 


set to 0.006. The right panel shows the distribution of the residual fitted by power function. 
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图 5 左 子 图 展示 了 2021 年 1 月 15 日 资料 (高 度 角 为 83~78”， 仪 器 星 等 为 11-13mag) 恒 星相 对 像素 距离 与 相对 测 光 
精度 之 间 的 关系 ， 其 中 剔除 阔 值 设置 为 0.005， 右 子 图 是 使 用 过 函 数 拟 合 后 残 差 随 相 对 像素 距离 之 间 的 分 布 情况 。 


Fig.5 The left panel shows the relationship between the relative distance and photometric precision of the observations on 


January 15, 2021 (the altitude is 83 ~ 78 deg and instrumental magnitude is 11-13mag), in which the elimination threshold is 


S set to 0.005. The right panel shows the distribution of the residual fitted by power function. 

© 从 上 述 结果 我 们 可 以 发 现在 同一 天 的 测 光 精度 来 看 ， 目 标的 高 度 角 越 大 ， 相 对 测 光 的 精度 越 高 。 

™ 在 天 顶 附近 时 ， 亮 于 13 等 的 恒星 测 光 精度 可 达 0.002 星 等 。 当 目标 高 度 角 较 小 时 ， 相 对 测 光 的 精度 

S 较 低 ， 且 存在 比较 明显 的 测 光 精度 回报 现象 。 从 图 3 中 我 们 可 以 看 到 ， 相 对 测 光 的 精度 随 相 对 像素 

69 距离 的 增加 而 降低 ， 测 光 精 度 从 0.003 星 等 降低 到 0.006 星 等 。 

N 3.3 精度 回报 现象 讨论 

SS 星光 经 过 不 同 质量 的 大 气 会 受到 不 同 的 消光 效应 。 因 为 相对 角 距 离 较 小 的 两 颗 星 经 过 大 气 的 质 

> 量 相似 ， 所 以 做 相对 测 光 时 存在 精度 回报 现象 。 当 目标 高 度 角 较 低 时 ， 大 气质 量 较 大 ， 相 对 测 光 的 

>< 精度 回报 现象 较为 明显 。 

= 在 地 面 观测 中 ， 通 常 我 们 会 选择 目标 在 天 项 方向 附近 观测 ， 以 减少 大 气 消光 的 影响 。 然 而 ， 有 

r= 些 目标 经 过 子午 线 时 仍然 处 于 较 低 的 高 度 角 ， 此 时 使 用 相近 的 参考 星 进行 相对 测 光 会 产生 较 大 的 收 

= 益 。 当 天 项 方向 的 大 气质 量 较 大 时 ， 可 能 也 会 存在 较为 明显 的 精度 回报 现象 ， 使 用 相近 的 恒星 作为 
参考 星 进行 相对 测 光 也 许 是 一 种 较 好 的 补救 方法 。 测 光 精 度 回报 现象 的 存在 引导 着 我 们 应 该 尽 可 能 


寻找 附近 的 孤立 亮 星 作为 参考 星 ， 当 然 我 们 也 需要 保证 参考 星 的 光度 是 稳定 的 。 
4 总 结 与 展望 


本 文 介 绍 了 现 有 的 测 光 工具 与 技术 ,提出 了 把 Maxim DL 软件 与 Python 语言 中 photutils 库 函 数 
相 结合 的 测 光 技术 。 通 过 该 技术 ， 我 们 能 够 方便 地 处 理由 于 星 像 半 高 全 宽 有 较 大 变化 时 导致 部 分 图 
像 目 标 不 能 被 检测 的 问题 ， 也 能 够 快速 地 进行 星 像 测 光 数 据 的 匹配 和 统计 。 通 过 分 析 中 国 科学 院 云 
南天 文 台 lm 光学 望远镜 在 不 同 高 度 角 拍摄 的 M35 星团 图 像 ， 本 文 发 现 了 在 目标 高 度 角 较 低 时 有 较 
为 明显 的 测 光 精度 回报 现象 。 在 后 续 的 工作 中 ， 我 们 将 结合 位 置 测量 和 光度 测量 的 精度 回报 现象 ， 
进一步 探究 天 体 中 光度 和 位 置 变 化 的 联系 。 
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Technology Research and Precision Premium in Relative Photometry Based on CCD Images 
Guo Bifeng'”, Peng Qingyu'*”, Shang Yijia'’, Qin Yonggui'”, Lin Furong'” 


(1 Department of Computer Science, Jinan University, Guangzhou, 510632, China) 

(2 Sino-French Joint Laboratory for Astrometry, Dynamics and Space Science, Jinan University, Guangzhou, 510632, China) 
Abstract: The changes of celestial bodies’ light with time can reflect their physical properties, and precise 
photometry technology helps us analyze and explore the changes of light better, which is of great significance to 
the study of the evolution of celestial bodies. This paper introduces the related photometry tools and technologies 
and suggests a method of photometry technology and its application by combining MaxIm DL with photutils. 
Besides, this paper puts forward a fast-matching solution of each star image’s photometric data in continuous 
frames. By analyzing the images of M35 cluster taken by the lm optical telescope of Yunnan Observatory, 
Chinese Academy of Sciences, we find there is a phenomenon of precision premium in relative photometry. 
Namely, two close stars have better measurement precision in relative photometry, and the significance of this 
phenomenon is related to the atmosphere mass. According to the phenomenon, this paper proposes some 


reference methods in high precision relative photometry. 


Key words: Relative Photometry; Precision Premium in Photometry; CCD Image Processing; photutils; MaxIm DL 


参考 文献 
[1] Devogéle M., TangaP., Bendjoya P.,et al., Shape and spin determination of Barbarian asteroids[J]. Astronomy & Astrophysics, 2017, 
607:A119. 
[2] Verbiscerl A.J., Porter S., Benecchi S.D., et al., Phase Curves from the Kuiper Belt: Photometric Properties of Distant Kuiper Belt 
Objects Observed by New Horizons[J]. The Astronomical Journal, 2019, 158(3):1-17. 
BI RHA, Sith, RÉE, 徐 小 林 . 平 谱 射 电 类 星体 3C 454.3 的 中 长 周期 光 变 特性 研究 下 . 天 文 研 究 与 技术 , 2020, 17(1): 1-7. 
WU Yuecheng, ZHANG Haojing, YU Lian, et al. The Medium and Long Period Light Variation Characteristics of FSRQ 3C 454.3[J]. 


Astronomical Research and Technology, 2020, 17(1): 1-7. 

[4] Landolt A.U., UBVRI PHOTOMETRIC STANDARD STARS AROUND THE CELESTIAL EQUATOR: UPDATES AND 
ADDITIONS[J]. The Astronomical Journal, 2009, 137:4186-4269. 

[5] Pascu. D., An appraisal of the USNO program of photographic astrometry of bright planetary satellites[C], Galatic and Solar System 
Optical Astrometry, 1994:304-311. 

[6] Morgado B., Assafin M., Vieira-Martins R., et al., Astrometry of mutual approximations between naturalsatellites. Application to the 
Galilean moons[J], Monthly Notices of the Royal Astronomical Society, 2016, 460(4):4086-4097. 

[7] LinF.R., Peng J. H., Zheng Z. J., et al., Characterization of the precision premium in astrometry[J]. Monthly Notices of the Royal 
Astronomical Society, 2019, 490(3):4382-4387. 

[8] Laher R. R., Gorjian V., Rebull L. M., et al., Aperture Photometry Tool[J], Publications of the Astronomical Society of the Pacific, 
2012, 124:737-763. 

[9] WEE, EE, 林 孕 荣 , Astrometrica 与 MaxIm DL 在 天 体 测量 的 应 用 与 技术 研究 [J/OL], 2020， 天 文 研究 与 技术 ， 
https://doi.org/10.14005/j.cnki.issn 1672-7673.20200701.001 


GUO Bifeng, PENG Qingyu, LIN Furong, Applications and Technology Research for Astrometrica and MaxIm DL in 
Astrometry[J/OL], 2020, Astronomical Research & Technology, https://doi.org/10.14005/j.cnk1.issn1672-7673.20200701.001 
[10] Da Costa G. S., Basic Photometry Techniques[C], Astronomical CCD observing and reduction techniques, 1992(23):90-104. 


